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Abstract 

Some rare decay processes are particularly sensitive to physics beyond the Standard 
Model (SM) because they have no SM tree contributions. We focus on one of these, 

— > 4>K S . Our study is in terms of the high scale effective theory, and high scale 
models for the underlying theory, while previous studies have been focusing on the low 
scale effective Lagrangian. We examine phenomenologically the high scale parameter 
space with full calculations, but largely report the results in terms of mass insertion 
techniques since they are then easily pictured. We also determine the ranges of different 
mass insertions that could produce large non-SM CP effects. Then we exhibit classes of 
high scale models that can or cannot provide large non-SM CP effects, thus demonstrating 
that data on — > <j)K s can probe both supersymmetry breaking and the underlying high 
scale theory and even make relatively direct contact with string-motivated models. We 
provide a novel and systematic technique to understand the relations between high and 
low scale parameters from RGE running. We include all constraints from other data, 
particularly b — > S7 and EDMs. 
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1 Introduction 



CP violation is a fascinating subject in particle physics. Complex parameters in the Lagrangian 
could give rise to numerous interesting observables in low energy experiments. Many, especially 
those associated with the third generation of quarks, are either not well measured or controver- 
sial. They represent new opportunities to discover new physics and uncover new flavor physics 
in the near future. If possible FCNC parameters are small, some presently unknown mecha- 
nism or symmetry has to be found to explain why they are small. Sizable new CP violations, 
as well as the KM phase itself, probe fundamental flavor physics, which is closely related to 
supersymmetry breaking and string theory. 

In the Standard Model, the KM phase 8km and the strong phase 6qcd are the only sources 
for CP violation. The bounds on the electric dipole moment (EDM) of the neutron imply that 
Oqcd < 10~ 9 , so we will ignore it for b physics. Therefore, effectively the only source of CP 
violation in the SM is the non-zero 8km- Many experiments have been analyzed to measure this 
phase and consistent results would increasingly establish the KM mechanism of CP violation. 
Before the measurement of the time dependent CP asymmetry in the process of Bd —> <f)Ks , 
all experiments indeed gave a consistent measure of 8km = 60° ± 14° pQ. On the other hand, 
many proposals for new physics beyond the Standard Model allow new sources of CP violation. 
Some of them are constrained by the experimental bounds from electric dipole moments and 
measurements in Kaon physics. Others could give rise to potentially interesting deviations. 

Low energy supersymmetry is the most compelling candidate for the new physics [2]. The 
most general soft supersymmetry breaking Lagrangian carries a significant number of phases 
which could be interesting new sources of the CP violation. Therefore, it is important to study 
the ways in which CP violation in the soft supersymmetry breaking Lagrangian could show up 
in the experiments. 

Since processes involving the superpartners (hence the new parameters) necessarily occur 
at loop level, the new CP violating effects could only be significant in the processes where SM 
contributions are also suppressed. Important examples of such processes include b — > 57 and 
B d -> <j)K s . 

It is interesting to notice that recently there is a potential discrepancy between the CP 
asymmetry measured in the process B d — > 4>K S (denoted by S^Ks) and the SM prediction: 

S+ K ~ S itK ~ sin 2(3 = 0.736 (1) 

The latest results from BaBar E] is 

S^ fl = 0.47±0.34i8;gg, (2) 

and from BELLE HE] is 

S+k 3 = -0-96 ± 0.50i£??. (3) 

Because this process is one of a few that are unusually sensitive to CP physics beyond the SM, 
it is worthwhile studying it in detail however the experimental situation is finally resolved. And 
of course if a deviation from the SM prediction is confirmed it is exceptionally important. 
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Since 5km is the only source of CP violation in the SM, if a deviation from the SM prediction 
for S^k is confirmed, it will be a clear sign of new physics beyond the SM and the new physics 
must be relevant to the weak scale. The superpartners in the loops must be at a mass scale 
that guarantees their production at the Tevatron and LHC. In the Minimal Supersymmetric 
Standard Model(MSSM), there are 42 new CP violation phases in addition to 5km and Qqcd- 
Although many of these phases are already constrained by various experiments there are 
still plenty of them which are not strongly constrained. It was shown in [TJ |H1 |H1 E3 EH > that by 
tuning some of these not strongly constrained CP violation phases at the weak scale, BELLE's 
result could be described without violating other experimental bounds. 

The most stringent contraints on CP violation arise from several electric dipole moment 
measurements. Since we are mostly interested in CP violation associated with the last genera- 
tions of quarks, the most relevant constraint will come from the mercury EDM which depends 
on the chromo-magnetic operator of the strange quark (see detailed discussion in Sec. EJ). From 
the point of view of a general low energy MSSM, this should not give a direct constraint on 
CP violation observed in flavor changing B decays since the EDM is only sensitive to flavor 
diagonal phases. However, we have to bear in mind that the soft Lagrangian is actually defined 
at some high scale at which supersymmetry is broken. Generically, we expect soft parameters 
at that scale to carry large phases that can give large deviations from the Standard Model in 
CP violating observables. The renormalization group running mixes the flavor- diagonal and 
flavor off-diagonal parameters and the mixing is enhanced by the large Logs in the conventional 
picture of gauge coupling unification. One of the main results of this paper is the study of this 
effect and its implications for the high energy allowed parameter space. 

With better understanding of the allowed high energy parameter space, we then study 
and constrain some models of flavor structures. Previous studies have focused on weak scale 
phenomenology. Here instead we emphasize the properties of the high scale theory and their 
relation to the data. We do both a phenomenological analysis of what high scale properties 
would be needed to explain a deviation from the SM, and an examination of what kinds of high 
scale underlying theories could or could not have the needed properties. 

This paper is organized as follows. We describes the most relevant CP violation observables 
and constraints in Section 2. A model independent study of the input high scale parameter 
space is performed in Section 3, utilizing both semi-analytic and numerical methods to study 
the RGE effects. With the knowledge obtained from the model independent study, we examine 
several classes of models of high scale flavor structure in Section 4. We present our conclusions 
in Section 5 and a detailed semi-analytic study of MSSM RGE effects using a new approach in 
the appendix, which provides a clear picture of the interplays among the flavor parameters. 

2 CP Violation Observables and Constraints 

In this section, we review the main experimental contraints and prospects of detecting CP 
violation beyond the Standard Model. We discuss the Mercury EDM and b — > 57 as they 
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usually provide the most useful constraints. We use Bd — > 4>Ks as the main example of an 
observable potential deviation. Other constraints jHj involving first or second families of quarks 
are less important for our purpose. We will not discuss them here, but still require our results 
to satisfy those experimental bounds. We emphasize that, although mass insertions generally 
provide a clearer picture of the flavor changing, we make use of the full squark mass matrix 
in our calculation. The necessity is evident in the discussions of the importance of various 
multiple-mass-insertions in the following. 

2.1 Mercury Electric Dipole Moment 

The EDM bound on Hg puts constraints on the CP violation phases of the soft terms. It 
requires the strange quark chromo EDM (CEDM) to satisfy ^21 H~3| 

|ed? | < 5.8 x 10~ 25 ecm, (4) 

assuming vanishing up and down quarks CEDM. However, there are significant theoretical 
hardronic uncertainties in extracting d% from Hg 199 1 . 

In the MSSM, although chargino exchange diagrams can give sizeable contributions to the 
strange quark CEDM [Tl], the contribution is usually dominated by the gluino exchange dia- 
grams 

e( f = ce—^lm(A d 2 ' 2 LR ) x L(x) + L <-> R, (5) 
An vnA 

with c = 0.91 |15j . L(x), where x = m 2 ~/m^ is a loop function (or its appropriate derivative). 
Notice it is different from c = 3.3 usually quoted in the literature [TSJ. This difference comes 
from the different definition of the chromo-dipole operator. The chromo-dipole operator used 
in obtaining eq.(EJ) is defined as | g s q~it a G a a^^qi which includes the strong coupling g s . With 
this definition, which is also used in a recent study ^7j, there should be no large scaling of this 
operator. 

A{j LR , {f = u -, d) is a generic mixing parameter between left-handed i-th generation and 
right-handed j-th generation of up- or down-type squarks. It could involve single or multiple 
mass insertion parameters which are directly defined from the soft parameters. Several examples 
are 

^22 LR = 7^(^)23(^)32, 7jj(^L/?)23(^Rij)23' 

&22 LR = 37(^)23(^)33(^)32, 37(^)23(^)33(^)32 (6) 

In the cases of double mass insertion, there is an extra factor of 1/2 coming from the Taylor 
expansion. Although flavor changing parameters do not contribute to the EDM at leading order, 

our study, we therefore allow the theoretical calculation to have a factor of 3 uncertainty. 



3 



they do enter at the next order through combinations with other flavor parameters. A particular 
triple mass insertion is studied in Ref f7j . Using Eq. we obtain for that combination 

edC ° = (-55) ^(^S&M^U, (7) 

where we have set gluino and average squark masses to be equal, rh g = m q , when evaluating 
the loop function. Notice also the extra 1/6 coming from the Taylor expansion. Combining 
this factor with the the scaling behavior, we found that the constraint of Ref. ^7] should be 
about a factor of 20 less restrictive than they report. 2 . Nevertheless, this combination will 
still provide some constraint on the high energy CP violating parameters. 



2.2 b S7 

b — > 57 is a process where the Standard Model tree level contribution is absent. Therefore, 
the SUSY contribution, which enters at one-loop order, could be comparable to the Standard 
Model processes. As a result, new flavor parameters in the soft Lagrangian will be significantly 
constrained by this process. We use the following bound in our calculations: 

2.0 x 10^ 4 < BR(6 -> s 7 ) < 4.5 x 10~ 4 (8) 

In the MSSM, the SUSY contribution to BR(6 — > 57) is usually assumed to be dominated 
by the chargino loop contribution to the 77 (or 0' 7 ) operator: 

New 2 m b tan(3 LL 

W 7 X 92 ^12^23 X L h _ w 

2 , 2 m t m 6 tan (3 , LR 

+ 92^ 2 A 33 V 22 X L h 



m 



+ 9 2 2^X 23 V n x L W . (9) 



where A 23 is either A^ L or A 2 . Vij are chargino mixing matrix elements. The three different 
lines correspond to Higgsino-Wino, Higgsino, and Wino loops, respectively. Typically, the 
Higgsino-Wino loop gives the dominant contribution. The others could be important in some 
circumstances as well. Eq. is schematic and designed to show the dependence on various 
flavor parameters. We have not shown explicitly the charged Higgs contribution as it does 
not depend on flavor changing soft parameters. There is a similar expression for the chirality 
flipped 0' 7 operator as well 

A™j AB s should again be understood as compounded parameters. In this paper, we are mainly 
interested in the flavor physics parameters associated with the down sector. However, as we 
will argue in section 4, some of them are related to the up-sector flavor parameter via 577(2) 

2 In a subsequent study |18j . the same group of authors changed their result after private communication 
from us. Their later results agreed with the results presented in this paper. 
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gauge symmetry and CKM mixings. Therefore, Br(b — > 57) will provide important constraints 
in those cases. On the other hand, the current experiments do not impose strong constraints 
on the CP asymmetry of b — > 57. So we will not discuss it here. 



2.3 B d (J)Ks 

The time-dependent CP asymmetry measured in the process of Bd — > <pKs or Bd — > J/ipKs is 
expressed in the following way: 

«,(t) = r(g(t)^/)- W ^/) 

= C/ cos AM Bd t + S f sin AM Bd t 
where / = 0i*Cs or J/ depending on which process we are studying. Defining as usual 

fill 




we have 

1-IA/l 2 2Im(X f ) 
1 + |A,| 2 ' f 1 + IA 



C f - 1 I IX 12' 6 / - 1 7 |\ 12 ( 12 ) 



In the Standard Model, we have 



^ = 1 = 1 < 13 > 



^ d V u V*Vcd 



This is invariant under redefinitions of phases. 

The SUSY contributions to A/ can also be written in a phase rotation invariant fashion. 
To simplify the formula, we assume SUSY contributions only modify the decay amplitude 
significantly, not the mixing part. We obtain 



v ,. = (14) 



where we have defined 



r ^ _ (A d 23 LR ) M 3 , 32 _ 632 (A% LR y M 3 * 

in order to distinguish the 0& g and Og 9 contributions. ^ is the ratio of the magnitudes of 
the SUSY and SM contributions. 6*23 and #32 are explicitly phase-rotation invariant. They are 
two of the 42 new CP violation phases in the MSSMj2U]. In other words, a non-zero #23 or #32 
characterizes CP violation beyond the KM mechanism. From eq. (H4j) . it's clear that when #23 
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and 632 are both zero, the time dependent CP asymmetry measured in the B factories should 
be 

S+k = sin 2(3 = 0.736 (16) 

and thus confirm the standard model KM mechanism of CP violation. Any deviation from 
this result would imply a new CP violation source other than 8km- in the framework of the 
MSSM, it would imply that some CP violation phase such as #23 or #32 is non-zero. Notice 
that although the SM CP violation phase f3 and the new phase(s) #23 (or #32) are independent 
parameters from the low energy effective theory point of view, it's still possible that these 
phases are correlated in fundamental flavor physics. Generically, the number of independent 
phases in the fundamental theory could be less than the number of phases carried by low energy 
parameters. 

In the Wolfenstein parameterization, Vt s and Vtb are real. In the MSSM, the gluino phase 
can be rotated to zero by using the U(1)r symmetry. Then #23 and 6*32 are just the phases 
of A^^ and A 32 R , respectively. In the following, we will use this parameterization, but one 
should keep in mind that when we say there is a non-vanishing phase of these two Mis, what 
we really mean is the phase of the combined quantity in eq.([15|). 

It was shown in jH] that if at the weak scale, 3 

A% LR \ W or A% LR \ W ~ O(10~ 2 ) x e * (17) 

where is a non-trivial phase, then S^k can deviate significantly from its SM value. We will 
follow Ref. jH] to calculate the hadronic matrix elements, using the BBNS method j2J- in 
this approach, the strong phases arise from four classes of diagrams, vertex corrections, pen- 
guins, hard scattering with spectator quarks and annihilation diagrams. It is important to 
properly account for the power corrections originated from the latter two classes of diagrams. 
These contributions should be subleading in the BBNS factorization, but they involve infrared 
divergent integrals. To regularize them, we follow BBNS and parameterize the integrals as 
A = (1 + pe 1 ^) log(ms/A^) with = 500 MeV. This parameterization introduces hadronic 
uncertainties into our calculation. However, as discussed in greater details in [S], the uncer- 
tainties affect the branching ratios more than the CP asymmetries, and we focus on the CP 
asymmetries in the present paper. The uncertainties also decrease quickly for heavier gluino 
masses if we assume the validity of BBNS factorization and choose moderate values of p. We 
henceforth set p = since our gluino will be heavy. In the present paper we also mainly 
compare CP asymmetries for different models, and we expect relative results to be even less 
sensitive to hadronic physics. 

In this paper, we would like to study high scale models which can give large non-SM CP 
violation in — > <pKs ■ We will use these as examples of classes of models which could give 
rise to interesting low energy CP violation beyond the SM. One interesting result is that some 

3 Hereafter, when necessary we will use (...)|w and (...)|a to denote the weak scale and GUT scale values, 
respectively. 



6 



classes of high scale theories do not allow one to obtain results such as eq. (|T7jl without violating 
other constraints. Thus the low scale data can rather directly probe high scale theories. 

3 Model Independent Study of High Energy Parameter 
Space 

3.1 RGE Running of High Energy FCNC Parameters 

Super symmetry is used to stabilize the large hierarchy between the GUT scale and weak scale. 
Therefore, the soft parameters at the input scale, which is assumed to be close to the GUT 
scale, are related to those in the low energy effective soft Lagrangian only after the inclusion of 
radiative corrections enhanced by large logarithms. For CP violation in FCNC processes it is 
crucial to take those RGE running effects into account. In this section, we briefly summarize 
the results which are very useful for qualitatively understanding the constraints on the high 
energy parameter space (see the Appendix for a detailed study). 

As explained in the Appendix, a systematic and novel way to study the RGE running of 
the flavor parameters is to use what we call the High Energy SuperCKM (HES) basis. In the 
HES basis, we are dealing with approximate physical parameters and enjoying the presence of 
several small parameters, such as A ~ 0.22. The qualitative result of that study is summarized 
in Table 



Parameter 


Universal Contribution 


Feeds into 




~ 7]\ 2 y 2 t ~ 0.01 




(<W) 23 


~ rfylvtX 2 < 10" 4 






~ ^vX 2 Vt < io- 4 


(S d LL )2S ~ 50(^)^(^)33 


(^1r)32 




(SUhs ~ 100(^)32(^)33 



Table 1: RGE Analysis of High Energy FCNC Parameters. 77 ~ \tsw — £gc/t|/167t 2 ~ 0.2 is the 
loop integration parameter. 

Notice that starting from a universal boundary condition at the high scale, the FCNC 
parameters still acquire non-zero values due to small CKM mixing effects enhanced by the 
RGE running. We call them universal contributions to the flavor parameters. Important 
effects of this type are summarized in the second column of Table ^ They should be regarded 
as important "model-independent" values of those parameters 4 . 

4 Of cause, one could start with a set of non-zero FCNC parameters in such a way that their initial values 
exactly cancel the RGE generated contribution. This would require a conspiracy between the high energy 
fundamental flavor physics and the radiative corrections associated with lower energy scales. We do not consider 
such a possibility. 
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It is also important to notice that different FCNC parameters generically mix through the 
RGE running. In particular, starting with one FCNC parameter, others will be generated 
through the mixing. We document some of the most important mixing effects in the third 
column of Table ^ One of the obvious features is the "decoupling" bbehavior of LL and RR 
flavor parameters. Being dimension two parameters, they would not enter the running of the 
trilinears. As discussed in the appendix, the mixings between LL and RR flavor parameters 
are also suppressed either by second generation quark masses or by small CKM mixings. 

Using the results of the RGE study, we now examine the connection of the high energy 
FCNC soft parameters and the low energy observables. We again remark that while we present 
results for simplicity i terms of mass insertions, we actually do full numerical analyses without 
using the mass insertion approximation. 

3.2 High Scale Parameter Space 

We focus on the three observables discussed in the previous section. In a generic flavor model, 
all of the flavor parameters could be non-zero at high scale. However, in order to illustrate the 
constraints effectively, we study the cases in which only one of those parameters is non-zero at 
a time. The scenario is then referred to by the sole non- vanishing input scale flavor parameter, 
e.g., {8ll)2?\k GU t scenario which implies that (5^)23 is the only non- vanishing mass insertion 
at input scale A GUT . We will identify the input scale as the GUT scale, where the gaugino 
masses are 

Mi = M 2 = M 3 = 300 GeV, (18) 

unless otherwise noted. The overall scale of the squark mass matrices m and trilinear terms 
are set to be 300 GeV as well. We will choose tan/5 = 15 at weak scale throughout this section. 

3.2.1 (<5ll)23|a gc/t 

(^ll)23 could contribute to B d — > 4>K S through the chromo-dipole operator 8g 

°s 9 = -|^™&^(1 + l5 )G^b (19) 

with Cs g oc (5ll)23(^l_r)33- If it had a large CP violating phase it could give rise to a large 
deviation from the SM value of the CP asymmetry of — > <pKs ■ 

Unlike many other cases, this CP violating parameter is not constrained by the EDM bound 
at the high scale. The reason is as follows. The lowest order contribution of (5^)23 to the 
mercury EDM is proportional to (Slr)32(^ll)23- Using the result of Table ^ we see that 
starting with only (6^)23 at the high scale would not generate (5^)32 through RGE running. 
Therefore, this type of contribution to the EDM is suppressed. The next order contribution 
is through the combination (5^)23(5^)33(5^)23 oc (5^)23^^. As indicated in Table [TJ the 
universal contributions to (5^)23 due to RGE runnings is small. Adding the fact that LL 
and RR FCNC soft parameters do not mix significantly under the RGE running, this type of 
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Figure 1: The correlation of S^k s (solid lines) and BR(6 — > 57) (dashed lines) is shown in 
this plot. The x- and y- axises are the real and imaginary parts of (5^)23 at the input scale, 
respectively. The CEDM of the s quark receives relatively small contributions in this scenario 
and imposes no constraint. Here we set M2 = M3 at the input scale. 

contribution is also suppressed. To summarize, the scenario in which (6^)23 is the only flavor 
off-diagonal complex soft parameter at the high scale is almost unconstrained by low energy 
EDM measurements. Notice that the EDM constraint is usually the most stringent bound on 
the amount of the CP violation in the theory. Therefore, we consider (6^)23 as a promising 
candidate for large CP violation which could be probed in rare B-decay processes, such as 
B d - (j>K s . 

(611)23 is constrained by b — > S7 because it enters both the gluino and chargino diagrams 5 . 
Our numerical studies of (5^)23 ar e shown in Fig. [T] and Fig. El 6 We see that it is possible to 
satisfy the 6 — » 57 bound and generate an interesting CP violating effect in —>■ <pKs at the 
same time. More importantly, while giving interesting low energy CP violation this scenario is 
not constrained by the current EDM bound. For the entire parameter space shown in Fig. El 
|edf I < 8.5 x 10" 26 ecm. 

5 In the SCKM base, The contribution from chargino and up-type squark diagram involves the factor 
VcKMi^v^LLVcKM |22| . which is precisely the left-handed down type squark mass matrix (tti^)ll- Hence, 
(^ll)23 a l so contributes to b — > s-f through chargino diagram. 

6 We used leading order running of Wilson coefficients for all the calculations presented in this article. 
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Figure 2: The correlation of S^Ks (solid lines) and BR(b — > 57) (dashed lines) is shown in this 
plot. Here in contrast to Fig. ^ we scan M 2 (in GeV) and \{5 d LL ) 23 \ with M 3 = 300GeV and 
Arg((5f jL )23) = tt/2 fixed at the input scale. 

3.2.2 ( ^RR ) 23 1 J^GUT 

(S RR )23 could contribute to S^Ks through the 0' 8g operator as C' 8g oc (5f R ) l 3 (S RR ) 23- Since the 
right-handed rotations are not constrained by the Standard Model CKM matrix, (S RR )23 is not 
in principle related to {5 R r)23- Therefore, {5 RR )23 is not strongly constrained by b — > 57 . 

On the other hand, the mercury EDM strongly constrains the allowed CP violation carried 
by (5^)23- The leading order contribution is the combination (5^) 23 (5^) 2V The universal 
RGE contribution to (5^)23 is suppressed. Combining this with the fact that RGE evolution 
would not mix (5 RR )2 3 with (5^)23, we conclude the that leading order contribution to the EDM 
would not strongly constrain (5^)23- The next order contribution is (^ll)23(^l^)33(^_r)23 °^ 
(Sf L )23C 8g . As indicated in Tabled! there is a universal RGE contribution to (5^)23 ~ 0.01. 
Hence, a large non-SM CP violation in S^,k s from C' Sg would almost generically imply a large 
contribution to the mercury EDM. Therefore, the prospect of getting a large CP violation 
effect from (S RR )23, such as proposed in Ref. [Hj, is necessarily constrained by the mercury EDM 
bound, as indicated in Fig El 
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Figure 3: This figure shows the correlation between S^Ks (solid lines) and the s-quark CEDM 
(dashed lines in the unit of 10~ 25 ecm). (Sj lR )23 nearly does not contribute to b — > sj. In view 
of the uncertainties associated with the s-quark CEDM, we allow a relaxation of a factor of 3 
of the current experimental bound. 



3.2.3 (Slr)23\a gut 

($lr)23 could contribute to S^Ks through the 0$ g operator as C$ g oc (Sf R )23- It could give rise 
to a large deviation from the Standard Model predictions. 

The leading order contribution to the mercury EDM comes from combination (5^)23(^x^)23 oc 
(5 RR )2 3 C$ g . It does not strongly constrain (5^)23 since RGE running would not induce (5^)23 
either from (Sf R )23 or from universal contributions. The next order contribution comes from the 
combination (5x^)23(^^)33(^^)32- The universal contribution to (5| K ) 32 is highly suppressed. 
Since (5^)32 (approximately) does not mix with (6^)23 in RGE running, the contribution to 
the EDM is again suppressed at this order. Therefore, although (5^)23 does feed strongly 
into other soft parameters such as (5^ L )23, CP violation in (5^)23 would not be very strongly 
constrained by the mercury EDM. 

b — > s-f provides interesting constraints on (6^)23- For a large class of models, (5f jR )23irit/mb ~ 
^lr)23- ($t R )23 contributes to b -> S7 through the combinations (5^)23(^^)33 and ( s lr) 23(^)23- 

Numerical study of the (5^ R ) 2 3 scenario is shown in Fig. El We also include a study of the 
scenario where {Sf R )23 is not related to (5^)23, which is set to zero at the high scale. The 
result is shown in Fig0] 
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Figure 4: The interplay between S^Ks (solid lines) and BR(6 — > 57) (dashed lines) is shown. 
We assume that (#2^)23 = is not related to (5^)23 at the input scale. As a consequence, both 
BR(6 — > S7) and the s-quark CEDM become less constraining as opposed to the more usual 
case in Fig.© where Af' L = A® ,L at the input scale with A{j = Y/jA{' L being the trilinear 
term (/ = U,D). 



3 - 2 - 4 ($Lr)32\agut 

(5^)32 will contribute to S^,k s through the 0' 8g operator as C' 8g oc (Sfj^^- Since this is the 
first order in the mass insertion, generically, (5^)32 could give rise to larger CP violation in 
C 8g than (5^)23- According to Table [TJ an input scale (5^)32 also feeds strongly into (5^ R ) 2 3 
through RGE running. Therefore, it could induce a larger contribution to C' 8g through (5^)23- 

Since (5^)32 is not generically related to (6^)32 at the high scale (see Section HJ), b — > S7 
does not seriously constrain this scenario. 

Since the CP violation entering S^k s is from the operator 0' 8g , it is again strongly con- 
strained by the {S'[ L )23C 8g contribution to the EDM. The constraints are similar to those 
derived in the (5^)23 scenario. The numerical result is presented in Fig. The conclusion is 
that the prospect of CP violation in (5^)32 is highly constrained by mercury EDM bound. 
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Figure 5: This figure shows the interplay among S^Ks (black solid lines), BR(6 — > 57) (blue 
dash lines) and the s-quark chromo-EDM (red dash-dotted lines, in the unit of 10~ 25 ecm). We 
assume Af' L = A®' L at the input scale, where Af. = Y/jA{' L is the trilinear term (/ = U,D). 



4 Models of High Energy Flavor Structure 

In this section, we study several different models of high energy flavor structure. Our primary 
interest is in which models can have large low energy CP violation beyond the Standard Model 
and which ones cannot, for generic reasons. Given its considerable current interest, we will 
focus on the — > <pKs CP asymmetry. Our analysis could be generalized to other low energy 
CP violating observables. 

4.1 General Estimates from Supergravity 

In this section, we give general estimates of various mass insertion parameters using the general 
supergravity Lagrangian. To begin with, we discuss a parameterization of the trilinears with 
rather general assumptions about the supergravity induced soft terms. 
In supergravity based models, the trilinear terms can be written as 

Mj = Y%A% A% = Y?A%. (20) 

Notice there is no matrix product in these two equations. It is in this parameterization that 
the absolute values of A^- D are at the order of 0(m 3 / 2 ). One can certainly write trilinears as 
A D = Y D ■ A D where a matrix product is used. In this form, it's not guaranteed that every 
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Figure 6: The correlation between S^Ks (solid lines) and the s-quark CEDM (dashed lines, in 
the unit of 10~ 25 ecm) is shown as a contour plot. Br(B — ► Xg^) imposes no constraint here, 
while the s-quark CEDM is very constraining. We relax the CEDM bound by a factor of 3 to 
remind the reader of its large uncertainties. 

element of A is of order Oim^i-i) and the structure of A D could be quite different from A D in 
eq. (|2()jl . Since we assume the MSSM Lagrangian arises from supergravity theory, we will use 
the parameterization in eq. (|2()jl . 

To get a good estimate of S d ' LR , a reparameterization of eq.()20|) will be useful. Take A D as 
an example. It can be written as 

A% = Aq + A® ,L + Af' R + A^SaSji + A»5 iX 8 j2 + A^S i2 S n + A%6 i2 6 j2 (21) 

where 
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- A D - 

^23 
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+ ^33 



(22) 



If Yukawas are hierarchical, terms with a prime on the RHS of eq. (|21j) have to multiply small 
elements of the Yukawa matrix to give the corresponding elements in the trilinears A D . Thus 
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under the assumption that all Aij are of order the gravitino mass and the Yukawas are hier- 
archical, these terms are suppressed by at least m s /m b ~ A 2 where A = 0.22 is the Cabbibo 
angle. Hence we can safely neglect them. Keeping this in mind, one has 

A% « Y%(AV + A?' L + A? R ) (23) 
~A% - Y?(A° + A?> L + Af' R ) 

In most models, A^' L = Af' L = Af since u l L and d l L are in the same SU(2) doublet. Under 
this assumption, the above equations reduce to 

35 = Y%(AV + A? + A? R ) (24) 
A% = YP(AZ + Af + A?> R ) 

The authors of j2Hl showed that trilinears of many SUSY breaking models follow this parame- 
terization even without the assumption of hierarchical Yukawa couplings. Thus eq.([24|) can be 
applied to a wide range of models. 

To relate the trilinear terms to observables, we need to rotate A to the SCKM base. To be 
precise, let's define the down-type squark mass insertion(MI) as 

5 a, LR = (KEA D K%%v d 

( m LL)ii ( m RR,)jj 

where is the vev of the down type higgs and matrices K D are the rotation matrices that 
diagonalize Yukawa matrices. Among these Sfj LR s, the most interesting ones are 52 3 LR and S 32 
since they are directly related to the B — > 4>K$ signal. (Notice that 8^ = 8- * so we only 
need to study 5^ .) By using trilinears as written in eq.(j2U), we obtain 

= ,~ D , \~D s x (26) 

(m b (A^ L - A°> L )(KE) 23 (KE*) 33 + m b (A?> L - ^ L )(^£) 21 (^£*) 31 

+m s (Af' R ~ ^ R ){K°) 23 {K°*) 33 + m s (A?' R ~ ^ R ){K°) 21 {K°*) 31 ) 

Using the fact m s < m 6 , (K£) 23 (Kf*) 33 < (Kf ) 21 (Kf *) 31 and the definitions in eq.(j22) we 
simplify the above formula to 

rd,LR „ m b (-A°> L )(KE) 23 (KE*) 33 _ m b (Ag-Ag)(KE) 23 (Kf*) 33 

°23 ' 2 ~~ 2 

m~ m~ 

Using the same approximation, one also gets 

,a,L R „ M-< R )(K°) 33 (K°*) 23 _ m b (Ag-Ag)(KE) 33 (K°*) 23 

°32 ' o ~~ 2 y z ° 

Eq. ()27l28|) give a quite good estimate of #23^ an d 3 32 LR at the high scale. They suggest that by 
increasing the splitting of trilinears between the 2nd and 3rd family and/or increasing specific 
elements in the mixing(i^ D ), these two Mis could be enhanced. 

Next, we make numerical estimates of individual mass insertion parameters. 
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4.1.1 5 23 LR 

The general form is given in eq . (|27|) . The left-handed rotation iff is constrained not only by 
unitarity, but also by the CKM matrix: 

Vg KM = (KY) 21 (KE*U + {K U L ) 22 {K°*) 32 + (KY) 23 (K£ *) 33 = A 2 (29) 

where A = 0.22 is the Cabibbo angle. If there are no large cancellations among the 3 terms 
in the middle part of the above formula, one expects that each of them is less than or around 
A 2 ~ 0.04. This suggests (iff ) 23 is much less than 1 and implies a large suppression on 5 23 LR . 
To avoid this suppression, one has to assume both of the last two terms in the above equation 
are large and they somehow cancel each other to get a small number A 2 . This could happen in 
models with democratic Yukawa couplings. If this cancellation indeed happens, then from the 
unitarity constraint: 

I(^f)l3| 2 + |(^f) 2 3| 2 +|(^f)33| 2 = l (30) 

the maximum mixing we can have is |(iff ) 23 (iff *) 33 | m 0.5. In gravity-mediated SUSY break- 
ing models, it's natural to assume that both A D and nig are of order 0(m 3 / 2 ). Thus one gets 
S 23 LR « 0.5 x m b /m 3/2 . Suppose m 3/2 » 200GeV, then 5 23 LR « 0(0.01). 

This estimation is made at the high scale. To make contact with observables, one must take 
the RGE running effect into account. The main RGE effect is to enhance the diagonal terms 
in the squark mass matrix due to the gluino contribution. The off-diagonal terms don't run 
muchj^Zj. Thus in eq.(|27|) and eq.(|2Sj). only the denominators are significantly affected by the 
RGE running. The RGE running of the diagonal squark masses is approximately 

ml\ w ~ Qm\ /2 + ml (31) 



1/2 

6(v / 3m 3/2 )^ 
19m 2 /2 



™3/2 



On the second line we assumed a dilaton dominated SUSY breaking scenario: m\i 2 = V3m 3 / 2 
and m = m 3 / 2 . In the above formula, m\i 2 and m on the RHS should take their high scale 
values. This formula shows the low scale value of S 23 LR will get a factor of 19 suppression 
from its high scale value in the dilaton dominated SUSY breaking scenario, or a factor of 7 if 
one assumes vn\i 2 ~ m « m 3 / 2 . Thus the natural value of 5 23 LR is too small to give a large 
deviation from the SM for the B^ — > 4>Ks process. 

To compensate the RGE suppression, one needs to increase the splitting between A® 3 and 
Af^. For example, if we take A 23 = —Af^ = 4m 3 / 2 , there will be a factor of 8 enhancement in 
the numerator of eq. (|2*7j) . This allows 5 23 \w to be of order 0.01 which could give rise to large 
CP violation in — > <f>Ks if we assume the cancellation happens in eq.()30|). 

From the discussion above, we see that in order to use a large S 23 LR generating CP asymmetry 
in the Bd — > (j)Ks process, the following conditions should be satisfied: 

• Large mixing in iff: | (iff ) 23 (iff )* 33 \ ~ 0.5 
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• Large splitting between A 23 and A 33 . 
As we discussed, to satisfy the first condition, a large fine-tuning may be required. 



4.1.2 S 32 LR . 



d,LR 



The basic formula for 8 32 LR is shown in eq . (|28|) . The numerical estimate is similar to the 5 23 
case and there are two conditions to be satisfied if one uses 5 32 \\ to generate a large non-SM 
CP violation: 

• Large splitting between A 32 and A 33 : (A 32 — A 33 ) ~ 81713/2 

• Large right-handed mixing \{K^)* 23 {K^) 33 \ ~ 0.5 . 



Unlike (Kg ) 23 , there is no CKM constraints on (K§] 



23- 



4.1.3 Double Mass Insertions from 5 23 LL and 5 23 RR . 

At the weak scale, both single MI and double MI can give significant contribution to B — > (j)K$ 
process. In the MI approximation, a double MI: 5 33 LR x 8 32 has an effect similar to a single 
5 32 LR MI, and S 33 R x 5 23 has an effect similar to S 23 R . 

The size of 5 23 can be estimated as follows. Assuming the right-handed down-type squark 
mass matrix is diagonal in the gauge eigenstates, we have 

#» = (32) 
m 



3/2 



fm 2 



l z -ml 2 ){K°) 32 {K°*) 33 
1 



m l/2 



Notice there is no fermion mass suppression and the only restriction on is unitarity. For 
S 23 LL , a similar estimate gives 

x d,LL tV D\ < v d*\ I A2 without cancellation 

*23 - ( £ )23(^ )33 ~ I Q 5 w . th canceUation (33) 

As in the single mass insertion case, the above estimates are at the high scale. The RGE will 
induce large suppressions at the weak scale. 

4.2 Abelian Flavor Symmetry Models 

The discussion in the previous section is quite generic, not relying on any flavor models. In this 
section, we study abelian flavor symmetry models and discuss whether they can give rise to a 
non-SM CP violation. 
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Abelian flavor symmetry models are interesting since they give a nice explanation for the 
observed hierarchical structure of fermion masses and the CKM matrix. In the supergravity 
framework, the soft SUSY breaking terms are related to the quantum numbers of the flavor 
svmmetries|28|. Therefore from a flavor symmetry model, we can calculate the parameters of 
the soft supersymmetry breaking Lagrangian and learn whether such models can give large 
deviations from low scale SM CP violation. 

Suppose there is a U(l)x flavor gauge symmetry at the unification scale. A field is a SM 
singlet but carries U(l)x charge: = —1. Gauge invariance requires that the superpotential 
take the form: 



W = ^YP9( qi + dj + h d ) 

ij 

+ y £Yg6(q i + u j + h u ) 



qi+dj+h d 



Mr 



QiDjH d 



(34) 



qi+Uj+h d 



Here M P denotes the Planck scale. q iy dj and ha are the U(l)x charges for Q iy Dj and H d , 
respectively. are some 0(1) numbers. 8(x) = 1 for x > and otherwise. To get a 
hierarchical Yukawa matrix, (</>) jMp should be a small number. We assume it is approximately 
A ~ 0.22. By choosing the U(l)x charges correctly, one can generate hierarchical Yukawa 
matrices such as: 

I A 8 A 5 A 3 \ / A 4 A 3 A 3 \ 



Y l 



A 7 
A r > 



A 4 
A 2 



A 2 
1 



Y D oc 



A 3 
A 1 



A 2 
1 



A 2 
1 



(35) 



This set of Yukawas could give correct fermion masses and a correct CKM matrix. Our discus- 
sion below doesn't depend on the detailed structure of the Yukawas. 

If the superpotential is coming from heterotic string theory, modular invariance conditions 
should be satisfied. To specify these conditions, we first write down the Kahler potential for 
the moduli fields: 

JT = - 5>g(T a + 73 (36) 

a 

Here we denoted both T-type and U-type moduli fields collectively by T a . For each matter field 
$, we denote the modular weights as n|, corresponding to T a . In the superpotential, YV(Y®) 
may also depends on T a and have modular weights: nfj^n^^). To keep the theory modular 
transformation invariant, the following conditions should be satisfied [28J: 

(ft + dj + h d )n% + n a Qi + n% J + n% + n a Djij + 1 = (37) 

In the following, we would like to argue that in this model, it's difficult to give a large 
beyond the SM contribution to the B d 4>Ks process. From the discussion of the previous 
section, we know that to get a large S^ R , (K^) 2 3 should be at 0(1). This requires 32 and 33 
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entries in Y D have similar magnitude, which implies that they have the same powers of A due 
to eq.flSU), i.e. 

q 3 + d 2 + h d = q 3 + d 3 + h d d 2 = d 3 . (38) 

In the flavor symmetry models, one usually assumes in eq . (|35|) are 0(1) and tries to explain 
the hierarchical structure by using different powers of <j)/Mp. Under this assumption, should 
be independent of moduli fields T a , or depend on T a in the same form. Therefore, all are 
equal. Then using d 2 = d 3 , we have 

nl 2 = n% (39) 

The two relations: d 2 = d 3 and = n^ 3 have important implications for the soft terms. 
To calculate them we assume a diagonal Kahler metric for the observable sector fields but 
allow non-universality of the diagonal elements since different fields can have different modular 
weights. The scalar masses and trilinear terms take the form[ 

2 2 
l ij — m 3/2l 



/9 (l + 4 + 3cos 2 ^^<e')^ (40) 



A% = A°Yg = m 3/2 (-^ S m6+(q l + d J + h d ))Y t f (41) 

In the above equations, 9 is the Goldstino angle and @« — 1- F° r the right-handed down 
type squarks, the second and third generations have the same U(l)p charge and same modular 
weights, so we have 

m 2 D =m* (42) 



Thus eq.(J2EJ) should give a good estimate for 5 32 and we have 



d 2 ~ "»d 3 - 
;hat 

d,LR 



Therefore, S 23 RR = 0. For the trilinears, notice that they satisfy the parameterization of eq. 



jALR „ m d {A» - A°){K°) 33 {K»*) 23 

d 32 -g (43) 

m d (d 2 -d 3 )m 3/2 (K°) 33 (K%*) 23 







So in the large right-hand down type quark mixing case (Y 3 ® m Y 3 ^), we have S 23 RR « S 32 LR ~ 0. 
By the same argument, one can show that in the large left-hand down quark mixing case: 
Y23 ~ 5*33 , we have S 23 RR ~ S 23 LR ~ 0. Therefore, in the large mixing cases, it's difficult to 
generate large non-SM CP effects. 

One can also try models with 0{\) suppressed mixing. For example, Y 3 ® : Y 3 ® = O(X). 
Then (K^) 23 « A and according to the estimate in section 14.1.21 \5 32 R \ will be much less than 
0.01. In this case, using the formula in section 14.1.31 we get 1^23^1 ~ A at the GUT scale. After 
running down to the weak scale and therefore including a factor of about 7 RGE suppression, 
we get 15^1 ~ ~ 0.03. This MI contributes to B d — > <fiK s via a double mass insertion, 
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which involves 5^. To estimate S^, let's assume rriq — n — 600 GeV and tan/3 = 50. Then 
we have 

^ m^tan^ ^ 0.24 (44) 
where is the running b-quark mass. Thus we can estimate the double mass insertion 

^23 R x $33 R ~ 0.03 x 0.24 « 0.007 (45) 

Remember for the double mass insertion case, there is an extra 1/2 suppression from the Taylor 
expansion. In addition, for this case, the 2nd derivative of the loop function is smaller than 
the 1st derivative for the region around m~/m~ = 1. Thus there is no large deviation from 
the SM for the CP asymmetry in the Bd — > <f)Ks process. Similar arguments apply to the 
Kjg : Y33 — 0{X) cases. Therefore, without fine-tuning some parameters, we conclude that in 
the 0(A) suppressed mixing cases, it's also difficult to obtain a large non-SM Bd — > <pKs result. 
In summary, for Abelian flavor symmetry models, under the assumptions: 

• The coefficients Yff and YV in eq.([34j) are all 0(1) and T-moduli independent 

• Kahler metrics are diagonal(and allowed to be non-universal) for the matter fields 

it is unlikely to get a large non-SM contribution to the CP asymmetry of Bd — > <pKs process. 

4.3 Family Dependent Kahlar Potential 

In the previous section we saw that if a U(l)p flavor symmetry model is assumed, it's difficult 
to give a large deviation from the SM for the Bd — > 4>Ks process. Therefore, in this section, 
we won't specify a particular flavor symmetry model. Instead, we take the Yukawas as given 
parameters. Hence there are no direct relations between the Yukawas and the soft terms. Then 
using a family dependent (but still diagonal) Kahler potential, we will show that it's possible 
to get a large SUSY contribution to the CP asymmetry in Bd — > <$>Ks ■ 

We first give a model which has large mixing between the 2nd and 3rd generation right- 
handed down type quarks. By splitting A^ and Af, we obtain a large ^32^^ and therefore a 
large contribution to Bd — > <pKs- Then we give a model which has large mixing between the 
2nd and 3rd generation left-handed quarks. By splitting A^ and A®, we have a large S^ 11 and 
it also gives a large contribution to Bd — > 4>Ks- In the large 8 3 2 case, the EDM bound puts 
strong constraints on how big the SUSY contribution to Bd — > <f)Ks can be, and in the large 
^23 LR case b — > 57 constrains it. 

It's also worth pointing out that in both models, the source of the CP violation phases 
reside in the Yukawas and all the SUSY breaking F terms(Fs and Ft) are real. If one takes 
the MSSM as an effective field theory, only the invariant phases, such as the SM CP violation 
phase 5km or SUSY CP violation phase in eq.()15|). are physical phases so that it doesn't matter 
whether we put the SUSY phases in the Yukawas or F terms. But from the underlying theory 
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point of view, Yukawas and F terms have different origins and represent different physics. Thus 
it's quite interesting that in our models, the only source of CP violation is in Y u and Yd and 
this may have important impact on string model building. 



4.3.1 Large 5^ LR case. 

In this model, we take the Yukawas as given parameters. They have the following form at the 
high scale: 



Y u = a u x diag{m u , m c , m t } 

Y d = a d x Vckm ■ diag{m d , m s , m b } ■ U 



(46) 



In these equations, a u and a>d are normalization factors depending on tan/5. Vckm and the 
quark masses should take their high scale values. U is a matrix which generate a right handed 
down-type quark mixing: 

/ 1 \ 



U 



sin uj 



V 



(47) 



/ 



cosu; 

— sin uj cos uj 

Notice if we don't have soft susy breaking terms, the phase <fi in U is not observable and can be 
rotated away by field redefinition. But with soft terms, certain combinations of <f> and phases 
in soft terms are physical observables and can't be rotated away. Also notice we assume the 
above Yukawas are for already canonically normalized fields. 

To calculate the soft terms, we assume a mixed dilaton/moduli susy breaking scenario and 
assume we are in the weakly coupled heterotic orbifold vacuum. The soft terms take the form 
as 



mi/2 
rrif 



m 3 / 2 sin I 



™3/2 1 1 



.4 



X 



3 cos 2 e^.Q 2 

-V3m 3/2 (sm6e^ s + 
1 + n? + n] + nl + (T ( 



cos 



-2-) 'a 



(48) 
(49) 

(50) 



a=l 

T*)d a log(Y ijk ) 



In these formulas, i,j,k denote different MSSM fields. 9 is the goldstino angle and @« = 1 
where a = 1, 2, 3 correspond to diagonal T-moduli fields associated with 3 compactified complex 
planes. For simplicity we assume only these moduli fields are relevant for soft term calculations. 
The 7s and 7 a are the phases for F$ and Fy a and we set them to be zero, so all the CP violation 
sources are in the Yukawas, n™ are the modular weights of a field % with respect to a— moduli. 
They are negative fractional numbers. For fields in the untwisted sector, the modular weights 
are = (—1,0,0) or (0, —1,0) or (0,0, —1) depending on which complex plane field is on. 
Modular weights for twisted fields are a little bit more complicated and actually we don't need 
to know them in this model. 
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Generally, the Yukawas depend on moduli fields so that the last term in eq.([HU|) is 
nonzero. Since we take the Yukawas as input numbers, there derivatives respect to T a can 
not be determined. Thus we have to assume their contributions to the trilinears are small and 
neglect them. Then the above formula for the trilinears satisfies the parameterization in eq. 

(j23). 

Using these equations we have 

^23 = A% — A3 



From the discussions in section |4~T1 we learned that to get large S^ 2 LR , ^23 should be large. Thus 
we need a small goldstino angle 9 which means a moduli dominated susy breaking scenario is 
preferred. 

In this model, we use the following parameters: 

1 7T 

m 3 / 2 = 420GeV sin# = -= u = - tan/3 = 40. (51) 

V7 4 

The relevant modular weights and Q a are shown in Tabled We assume the other matter fields 



Table 2: Qj and modular weights for some particles. 
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D 3 , 



have family independent modular weights. Therefore, all other fields have family independent 
trilinear and scalar masses. Notice n^ 2 and n^ 3 are different which means they have different 
Kahler potentials 

K ~ t 2 + t 2 * + t 3 + t 3 * + "' ^ 

This difference will only show up in the soft terms if cos 9 ^ 0. 

Taking these parameters as high scale input, we use RGEs to run the soft Lagrangian to 
the weak scale and calculate observables such as S^Ksi the strange quark CEDM, BR{b — > sj), 
the higgs mass, etc. We scan the phase (ft in eq.([47|) from to 2n. As explained before, for the 
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large 5 32 case, it's the EDM bound that is most difficult to satisfy. We show the correlation 
between predictions of S^Ks an d edf in figured The upper bound on \ed% | is 5.8 x 10~ 25 ecm. 
The theoretical prediction shown in this figure has a large theoretical uncertainty. If we allow 
factor 3 theoretical uncertainty, from this figure we find the smallest S^Ks * s about -0.3. We 
checked that for this model, other experimental constraints such as BR(b — > 57) and the higgs 
mass bound are satisfied. 

<,d,LR 

Large o^ 7 scenario 




'<|>Ks 



Figure 7: CEDM vs. S^k s m the large 5% 2 LR scenario. We scan the phase in eq. (}4Tj) from 
to 2ir. Two horizontal lines are the experimental bounds on the s-quark CEDM. If we use the 
exact CEDM bound, the smallest S^Ks we can S e t m this model is around 0.14. If we allow 
a factor 3 theoretical uncertainty on the CEDM prediction, the smallest SrpKs in this model is 
around -0.33. 



4.3.2 Large case. 

To make this MI large, we need a large mixing of the 2nd and 3rd generation left-handed quarks. 
Thus we use the following Yukawas: 

Y u = a u xU ■ diag{m u , m c , m t } 

Y d = a d xU ■ V C km ■ diag{m d , m s , m b } 

The parameters are: 

m 3/2 = 360GeV sin# = —= uo = \ tan/5 = 24 (53) 

v 7 7 

The a and modular weights are shown in Table El Due to the different modular weights for Q2 
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Table 3: Oj and modular weights for some particles 
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and Q 3 , the trilinears A® and are split. (Notice A Q = A U > L = A U > R .) From eq.(|27jl. we see 
that 52 3 LR will be large. We assume modular weights for other particles are family independent. 
We scan the phase <p in the U matrix. As explained before, in the large case, the SUSY 
contribution to S^,k s is mainly constrained by BR{b — > 57). We show the correlations between 
them in figure |H1 From this figure, we see that the smallest S^Ks we can get without violating 
the BR(b — »■ sj) bound is about 0.13. 

c d,LR 

Large o ? , scenario 




y — ' — ' — < — ' — 1 — ' — ' — ' — ' — 1 — ' — ' — ' — ' — l_ 
-\ -0.5 0.5 

^(|>Ks 



Figure 8: BRib — > 57) vs. S^Ks i n the large S 23 scenario. Two horizontal lines are the bounds 
on BR{b — > sj). The smallest S^,k s m this model without violating the b — > 57 constraint is 
around 0.13. 
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5 Conclusions 



Fundamental flavor physics is expected to manifest itself in the flavor structure of the soft super- 
symmetry Lagrangian at the supersymmetry breaking scale. In principle, there could be large 
CP violating parameters in the soft Lagrangian. They could give rise to large deviations from 
Standard Model prediction of the CP violating observables. One example is the recent potential 
deviation in S^k s of Bd — > <pKs ■ On the other hand, CP violation in the soft Lagrangian at 
high scale is constrained by low energy observables such as EDM and b — > 57. In this paper, 
we give a model independent analysis of the constraints on the high energy parameter space. 
Several interesting scenarios are identified which both satisfy the experimental constraints and 
give rise to large CP violation beyond the Standard Model. With the help of these results, we 
further investigated classes of high energy flavor structure in the soft parameters. We found 
that large classes of such models could not produce large deviation from SM in processes such 
as Bd — > <pKs and satisfy the constraints at the same time. Finally, we presented several 
scenarios where such a goal can be achieved. Thus such data could point toward some classes 
of high scale theories as favored. Further investigation along those directions are clearly very 
interesting and important especially if BELLE result on S^Ks °f Bd ~~ ^ 'P^s is confirmed. 
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A RGE study of flavor parameters 

In this section, we analyze the RGE running of FCNC soft parameters. The purpose of this 
study is to establish a systematic way to understand semi-analytically the running of flavor 
off-diagonal parameters in the MSSM. It provides an understanding complementary to that of 
the precise numerical study (which is presented in Section |3J) of the connection between high 
and low scale flavor parameters. The techniques of this appendix will be generally useful in 
future analysis. 

In general, the running of a specific parameter in the SUSY Lagrangian is basis-dependent 
since the RGEs depends on Yukawa couplings. The usual practice is to write the RGE equations 
in the gauge eigenstate basis (so that the gauge interactions are flavor diagonal). However, in 
this basis it is not possible to give a systematic and model-independent study of the running 
of the FCNC parameters due to the uncertainties of the Yukawa couplings. Therefore, it is 
more useful to study the running in the so called superCKM basis where there is a more direct 
connection of the low energy observables and high energy parameters. 

It is neither convenient, nor necessary, to diagonalize the Yukawa couplings at each stage 
of the running. In order to get a systematic estimate of the effect of the running, it is enough 
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to run the RGE in the superCKM basis defined by high energy input Yukawa couplings. More 
specifically, we diagonalize the Yukawa matrices at the high energy input scale. This diagonal- 
ization defines a high energy superCKM (HES) basis. We then run the RGE equations in this 
basis. Now, in this basis, the RGE equations are different from those written for the gauge 
eigenstates. However, it can be shown that the only new parameters showing up the RGEs 
are the CKM matrix (Vokm) obtained by diagonalizing the Yukawa matrices at the high scale. 
Due to the facts (which we will justify) that the CKM matrix has only small running [3U], and 
the diagonal Yukawa couplings in this basis are proportional to quark masses, it is thus possible 
to obtain reliable estimates of the running. 

We begin with Yukawa couplings. In the gauge eigenstate basis, the RGEs of Yukawa 
couplings are 7 

^ = ^[SYY^ + Y^ + STr^ 

^ = ^[^^SY^Y^ (54) 

where we have suppressed all leptonic Yukawa couplings. We then rotate to the HES basis by 
applying the transformation Y — > V£YV^ on both sides of Eq. |!3J We then obtain, in the HES 
basis (we still denote the Yukawa coupling in this basis by Y) 

it = + W/l^M^ + 3Tr(Fii;)i; - + 6g 2 2 + ^)FJ, 

dY 1 16 7 

~dt = 16^ ^CK M yuYlV^ KM Y d + 3Y d Y^Y u + 3Tr(Y^Y d )Y d - (-g 2 3 + 3g\ + -gl)Y d }{55) 

where the high energy CKM matrix Vqkm is defined by V^V d . In order to get an approximate 
solution to these RGE equations, we make the following assumptions 

1. The Yukawa matrices stay approximately diagonal in the running. Therefore, approx- 
imately, individual diagonal entries run independently and are proportional to quark 
masses. 

2. The off-diagonal running terms, which are proportional to Vq KM , can be treated as a 
perturbation. The main effect of such a perturbation is to generate off diagonal Yukawa 
couplings. 

3. The CKM matrix does not run very much (or the running effect is sub leading). Therefore, 
we have a systematic expansion (in terms of A) of the RGE effects. 

From these assumptions, we can solve for the additional mixing generated by the running 
(the RGE running of the diagonal terms is well known and dominated by SU(3) gauge coupling 

7 We are using the convention of Ref . ■ 
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and third generation Yukawas) by substituting diagonal Yukawa couplings in to Eq. and 
integrating approximately. 

As a result, we have the following approximate solutions 



(56) 
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(57) 



where rj ~ \tEW ~ ^gc/t|/167t 2 ~ 0.2. Since the diagonal Yukawa couplings are proportional 
to the quark masses, the most significant modification of the flavor mixing generated by RGE 



running are SV^ 
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0.1 A 3 and 5V 2 L 3 ~ 0.1A 2 . We see that the 13 and 23 elements of the CKM 
matrix, or a hierarchical V L proportional to the CKM matrix, could run about 10 percent. 
This is consistent with the assumption we made in solving the RGE equations. 

Next, we look at the running of the trilinears. In the HES basis, the RGE equation for the 
trilinear A d is 



dA a 
~dt 



1 



167P 



{4A d Y^Y d + 5Y}Y d A d + V^ WM Y u Y},V* KM A d + 2V^,,A,Y^ 
16 os „ .16 



'CKM 1 u 1 u v CKM d 



' 'CKM^u 1 u v CKmY<1 



+ A d (3Tr[Y}Y d ] - ^g 2 3 ) + Y d {f&M a + QTr [YjA d ])} 



(5? 



where we suppressed subleading terms (such as terms proportional to electroweak gauge cou- 
plings). First, we observe that the running of the diagonal terms of the trilinears is almost 
always dominated by the term proportion to the gluino mass (with the possible exception of 
a large 3rd generation diagonal trilinear coupling). Intuitively, the running of the off-diagonal 
terms are almost proportional to themselves. Therefore, their running should not be very sig- 
nificant. To gain an approximate understanding of the running, we could expand the RGE 
equations, as we have done in the case of the Yukawa couplings, in terms of small (off-diagonal) 
parameters such as A ~ 0.22. We write the RGE equation as 



167T 



, dA a 
~dt 



A d 



(59) 



where the righ-hand-side A d is a 3 x 3 matrix. In terms of small parameters, the flavor off- 
diagonal entries of the last two generations are 



(Adhs ~ -(y 2 + y t y b )\ 2 (A d ) 33 + (y t y b + 4y 2 b )(A d ) 23 + \y t y b (A d ) 13 
(A d ) 32 ~ 4Ay b {A d ) 33 + 5y 2 b {A d ) 32 , 



(60) 



where A ~ r]y b y^\ 2 , coming from term A d Y d Y\ in the RGE. The result of RGE running of 
off diagonal entries of trilinear couplings is approximately A d = r\A d . Of course, it is always 
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understood that the entries in matrix Ai should be taking some appropriate intermediate values. 
Useful estimates can be obtained from those expressions. For example, we can derive that even 
without the presence of off-diagonal terms in the trilinears, by RGE running, we will have 8 9 

fed \ v d v^ 2 yKA* d ) 3 3 ln _ 4 (au 

{(>lr)23 ~ — 10 , (61) 

and 



rrig nig 



d A7]Ay b v d (A* d ) 33 5 
{hR>32 —2 ~ 10 (62) 

m q 



Now we turn to consider the running of the soft masses. First, we consider the running of 
the right-handed down-type soft mass parameters. The RGEs, again in the HES basis, are 



! ^ = 16n 2 [2YjY d ml + 2miY+Y A + &lm\Y d + ^m\Y\Y d + AA\A d - 



We could write the RHS as a 3 x 3 matrix Aifj and expand it in terms of small parameters. 
For the last two generations, we have approximately 

(M 3 ) 23 ~ 2Ay h [(m & f m + (m 3 ) 2 22 + 2(m Q )l 3 } + 4(i d )* 2 (i d ) 33 . (64) 

Some important results can be derived from Eq. H3J First, if we begin at the input scale with 
a non-zero (Ad) 32, we would induce a right-handed mixing term through RGE running 

(4 R h3 - v A[AdM 2 AdU ~ A V (8i R )l 2 { -^%^ - 100(5^)^(^33- (65) 

Notice that a A 23 entry, on the other hand, does not generate a large RR mixing. This fact 
will be important in the search for viable high scale models. 

On the other hand, if there are only diagonal terms in the soft masses, the RGE running 
could generate an off-diagonal mixing 

(8 d RR ) 23 ~ 2r,Ay b [(m 3 )l 3 + (m D ) 2 22 + 2^)1,]/™? < 10" 4 (66) 

which is quite suppressed (comparing with the LL case studied below ). 

We also note that the mixing between (8^)23 and (5^)23 is highly suppressed by second 
generation quark masses and/or higher power of CKM mixing A. 

8 Notice that in our notation [2] , LR part of the squark mass matrix corresponds to A* 
9 In our estimation, we typically take the trilincar part to be the dominant part in the LR sector of the 
squark mass matrix. This assumption could be violated for the diagonal elements of LR sector, especially the 
down-sector 33 element, in the very large tan/3 and \i regime of the parameter space. In that case, the estimates 
proportional (or mass insertions proportional to (5f R )33 would be enhanced.) 
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Finally, the RGE running of the soft masses of the left-handed squarks, in the HES basis, 

are 

dm?- \ 

+ 2mH u VQ KM Y u Y u V£ KM + 2VQ KM A u A\y CKM 

+ Y d Y}ml d + m\Y d Y\ + Y d m\y\ + 2Y d m%Yj + 2A d A\ - -g 2 3 \M 3 \ 2 }, (67) 

which give 

(M Q ) 2 3 ~ -y 2 \ 2 [(m 2 Q ) 33 + (m|) 22 + 2(772^33 + 2m 2 H J 
+ Ay 6 [(m|)33 + 2(7711)33 + 2m 2 Hd ] + y 2 (m|) 23 
+ 2(A u ) 23 (A u y 33 + 2(A d ) 23 (A d y 33 + A(i u ) 13 (i u )* 3 . (68) 

Therefore, in the absence of any off-diagonal terms, the RGE running will generate a 

fe)23 ~ Wliimlhs + (m|) 22 + 2(m|) 33 + 2m%J/m% ~ 4 m 2 = 0.01. (69) 

Notice also, although a large A 2 3 does not generate a large RR mixing, it will generate a sizable 
LL mixing via RGE running. 

Similarly, the mixing between {5j jL ) 23 and (<5^r) 23 is highly suppressed as well. 
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